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ABSTRACT 



An experimental apparatus which will allow the study of 
the influence of ground effects on a heated cylinder placed 
in a crossflow of air was designed and constructed. 

The effects of a wall on the heat transfer character- 
istics of a heated cylinder were investigated at a Reynolds 
number of 153,000 for cylinder to wall spacings ranging from 
d/r = 0.0 to d/r = 5.33 (maximum), where d = gap width, and 
r = cylinder radius. Pressure coefficients were obtained at 
five degree intervals around the cylinder for each experiment. 

Experimental results obtained for d/r = 0.0 agreed well 
with previously reported investigations of a cylinder 
attached to a wall. Results obtained for d/r = 5.33 agreed 
well with previous work done with free cylinders. Experiments 
conducted at intermediate values of d/r resulted in minimum 
and maximum average Nusselt numbers, suggesting that gap 
spacing may exert a strong influence on the behavior of the 
flow field. 
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I. INTRODUCTION 



The principal objective of this work has been to design, 
construct and test an experimental apparatus which will allow 
the study of the influence of ground effects on a heated 
cylinder placed in a crossflow of air. Previous work of this 
nature has concentrated primarily on free cylinders. 

Recently, McComas (1) investigated the pressure and heat 
transfer characteristics of a heated cylinder attached to a 
wall. McComas determined the local and average Nusselt 
number, and the pressure coefficient, on the surface of the 
cylinder for Reynolds numbers varying from 53,000 to 153,000. 
This work is a follow-on study of the project initiated by 
McComas, differing in that the apparatus has been designed 
and constructed such that the plane surface can be varied 
from the attached position to approximately 3 cylinder 
diameters away. Due to the limitations of the wind tunnel 
used by McComas, (described in detail in (1)), hotwire 
anemometry and flow visualization techniques were not able 
to be used in his investigations. The present experi mental 
apparatus is located in the Mechanical Engineering Department 
wind tunnel in Building 500 on the Naval Postgraduate School 
campus. This wind tunnel, described in detail in section III, 
has a speed range of 10 to 175 miles per hour, and is equipped 
with a traversing mechanism and hotwire anemometer which can 
be used to obtain velocity information throughout the test 
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section. The tunnel also has a capability for flow visual- 
ization using smoke filaments. 

Experiments were conducted at a free stream Reynolds 
number of 153,000 for cylinder toplate spacings of d/r = 0.0, 
0.25, 0.5, 1.0, 2.0, 4.0, 5.33 (maximum separation) where 
d = gap width and r = culinder radius (Figure 1). The flow 
field in the vicinity of the cylinder was probed using the 
hot wire anemometer. The velocity distribution suggested 
that flow blockage was a factor in the experiment, and that 
the Reynolds number should be based on a corrected free 
stream velocity. The correction produced an effective 
Reynolds number of 161,000. Using this corrected Reynolds 
number, the data taken at d/r = 0.0 agrees with the results 
of McComas , and the data at d/r = 5.33 agrees with the work 
of Field (2). 

As was the case with McComas, experiments with the 
cylinder attached to the plate showed a critical flow pattern 
to exist over the top half of the cylinder. As the cylinder 
to plate spacing was increased, this flow pattern changed 
from critical to subcritical ever a narrow range of plate 
spacing (at about d/r = 0.09). The flow over the bottom 
half of the cylinder exhibited critical flow as soon as a 
space existed. This flow also changed from critical to 
subcritical, but the change was less pronounced and it was at 
a larger cylinder to plate spacing (at about d/r = 1.2). 
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Figure 1 . 



Schematic diagram showing determination of d/r. 



Another phenomenon which was seen to occur with increasing 
d/r was the gradual shift of the stagnation and laminar 
separation points over a range of 15 to 20 degrees. 

In calculating and plotting the average Nusselt number 
as a function of gap spacing, an interesting and unexpected 
trend was revealed. It was found that the minimum average 
Nusselt number did not occur when the plate was attached 
to the cylinder, but at a d/r = 0.25. Further, this curve 
also indicated a maximum point at d/r = 2.0. Although the 
current data are incomplete, (further investigation is 
needed over a range of Reynolds numbers) the Nu versus d/r 
curve at Re = 153,000 suggests that gap spacing may exert a 
strong influence on the behavior of the flow field. Future 
studies employing both flow visualization as well as heat 
transfer measurements would aid in establishing a plausible 
explanation for the trends observed in the curve of average 
Nusselt numbers. 
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II. BACKGROUND 



Although a great deal of investigation has been devoted 
to the flow and heat transfer characteristics of cylinders 
far from boundaries, see Meyer (5) and Field (2) for a list 
of references, virtually no references have been found 
relating to the heat transfer and flow characteristics of a 
cylinder placed in close proximity to a plane surface. The 
practical applications of such information would extend from 
aiding in conducting thermal analysis of container stored 
ordnance, rockets and missiles to the determination of the 
optimum size and spacing for heat exchanger tubes. 

Although McComas (1) experimentally determined local 
and average heat transfer coefficients on the surface of a 
uniformly heated cylinder attached to a wall, McComas found 
no references to either theoretical or experimental work on 
this subject. McComas proposed a theory, based on theories 
that had been postulated for forward and rear facing steps, 
that quantitatively explained the trends in his heat transfer 
data. Other than McComas' work, no further references have 
been found dealing with attached cyl i nders . 

Since the fluid flow field and thermal field are 
intimately coupled in a forced convection process, an under- 
standing of the flow field is mandatory if one is to under- 
stand the heat transfer process. A thorough discussion of 
the fluid flow and heat transfer characteristics of a free 



cylinder can be found in (3), and a proposed but as yet 
unverified description of the flow around an attached 
cylinder can be found in (1). 
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III. EXPERIMENTAL APPARATUS 



A. WIND TUNNEL 

All tests were conducted in a subsonic, open circuit, 
wind tunnel manufactured to Navy specifications by the 
Kenny Engineering Corp., of Monrovia, California (Figure 2). 
The overall length of the tunnel was 40 feet, with an 
entrance bell 5 feet, 5 inches wide by 7 feet, 5 inches high. 
The test section was 20 inches wide by 28 inches high at 
the entrance, expanding to 31 inches high to allow for 
boundary layer growth. These dimensions resulted in a 
contraction ratio, based on area, of 10:1. The tunnel was 
powered by a 75 HP, 1750 rpm variable speed motor driving a 
45 inch diameter adjustable pitch vaneaxial fan capable of 
supplying air at 10,000 to 70,000 cubic feet per minute. The 
tunnel was equipped with six graded mesh screens to control 
turbulence and a 24 inch manometer tube tapped into the 
stilling chamber and test section for velocity measurement. 
The tunnel was designed to operate at velocities ranging 
from 10 mph (14.67 ft/sec) to 175 mph (256.67 ft/sec). 
Turbulence intensity was taken using a hot wire anemometer 
and found to be approximately 0.2% at a free steam velocity 
of 98 ft/sec (see Appendix D). 
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Figure 2. Schematic diagram of the wind tunnel. 



B. HEAT TRANSFER EXPERIMENT 

The test cylinder was essentially three cylinders 
attached end to end (see Figure 3). The two end sections 
were hollow softwood cylinders supported by a hardwood 
brace, curved to match the cylinder contour. The central 
cylinder was a sheet of resistive paper (Temsheet) fitted 
into stepped portions of the end cylinders, and attached to 
the end cylinders and hardwood brace by double sided adhesive 
tape. The three cylinders then formed one cylinder 20 inches 
long, to correspond to the test section width, and having a 
diameter of 3 inches, chosen to allow maximum resolution of 
the liquid crystals without excessive blockage within the 
test section. The assembled test cylinder was tightly 
packed with glass wool to minimize heat transfer into the 
cylinder, and thereby reduce the time needed to reach steady 
state after a voltage was applied to the paper. 

The resistive paper used was a commercially available, 
carbon impregnated paper, Temsheet , with a nominal thickness 
of 0.039 inches. The width was sized to conform to the 
stepped circumference of the end cylinders (9.422 inches) and 
the length was 12 inches. The circumference of the resistive 
paper cylinder was marked in 5 degree increments, in two 
locations 2 inches apart. This was the test section used in 
all heat transfer experiments. The entire outer surface of 
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Figure 3. Photograph of the Temsheet cylinder. 



the resistive paper was covered with three coats of NCR S-43 
liquid crystals having event temperatures of: 

red: 109. 8°F 

green : 1 1 1 . 0° F 

blue: 112. 3°F (Field (2)) 

Power was applied to the paper by means of wires lead 

inside the cylinder through one end. The wires were attached 

to aluminum foil tape by means of alligator clips, and the 

edges of the tape were painted with silver paint to assure 

good electrical contact with the paper. The nominal distance 

between the tape/paint was 4.8125 inches, resulting in a 

2 

surface area for power dissipation of 45.34 in . 

The power source used was a Lambda regulated power 
supply, model LK 345A FM. The voltage applied by the Lambda 
was measured by a Fluke 8100A digital multimeter. This 
multimeter was also used to measure the resistance of the 
resistive paper. Air stream temperature was sensed by a 
copper-constantan thermocouple installed in the stilling 
chamber. The thermocouple output was read on a Leeds & 
Northrup millivolt potentiometer. 

The "wall" used in these experiments was manufactured 
from 1/2 inch plexiglass. This platewas46 inches long, the 
maximum practical length, to minimize the influence of the 
trailing edge on the flow patterns around the cylinder. The 
width of the plate was 20 inches, to conform to the test 
section width, and it was beveled on the under side of the 
leading edge at an angle of approximately 12 degrees. 
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The plate and cylinder were placed in the tunnel 
horizontally, and in such a way that the leading edge of the 
plate extended four cylinder diameters (12 inches) ahead of 
the cylinder centerline. The cylinder was attached to the 
test section walls by portions of the end cylinders pro- 
truding through holes cut for that purpose. The plate was 
capable of being moved vertically by means of a threaded 
rod extending through the tunnel floor, and a rear brace 
which locked in position by means of a split nut (Figure 4). 
The plate was capable of being located at any position 
relative to the cylinder, from attached, to a maximum 
possible distance of 8 inches away (Figure 1). 

C. PRESSURE COEFFICIENT EXPERIMENT 

The experiments conducted to find the variation of 
pressure around the ci rcumference of the cylinder used the 
previously described plexiglass plate as the wall and a 
3 inch diameter aluminum cylinder fitted with four pressure 
taps 90 degrees apart (Figure 5). The ci rcumference of the 
cylinder was scribed in 5 degree increments. Plastic tubing 
led from the taps through one end of the cylinder to a 
manometer bank. 

The manometer bank consisted of several tubes that were 
inclined at an angle of 30 degrees from the vertical. All 
four cylinder pressure taps as well as the stilling chamber 
and test section pressures were read from this manometer 
bank . 
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Figure 4. Photograph of the plate and cylinder assembly. 
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Figure 5. Photograph of the cylinder used to collect pressure data. 





IV. EXPERIMENTAL PROCEDURE 



The present work was conducted utilizing different test 
facilities than those used in previous investigations. The 
wind tunnel was smaller, and of the open circuit variety. 

Due to experimental considerations, the cylinder size and 
orientation were different, and the problem of varying the 
relative separation of the plate and cylinder was unique to 
this work. In order to assure that meaningful and accurate 
data would be obtained, it was necessary to compare initial 
results with previous work. Since both Field and McComas 
performed experiments at a Reynolds Number of 153,000, it was 
decided to check the present apparatus at that Reynolds 
Number. For this reason, an air stream speed of 98 ft/sec 
was used for comparison experiments. Because of time 
limitations, all of the data for the present work were taken 
at 98 ft/sec, corresponding to a Reynolds Number of 153,000, 
where Reynolds number is defined as: 



Re 



v 



Re = Reynolds number 



U 



free stream velocity at test section 
entrance, ft/sec 



co 



d 



cylinder diameter, ft 

kinematic viscosity of air, ft /sec 



v 
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A. HEAT TRANSFER EXPERIMENT 

The same experimental procedure used by Field and 
McComas was used in the present experiments. For complete- 
ness, a description of the procedures will be repeated here. 

Prior to the start of the experiment a sufficient 
voltage was applied to the cylinder to cause the liquid 
crystals to reach their event temperature and change color. 
This voltage was usually from 12 to 15 volts. After about 
a half hour, time to allow the cylinder and glass wool 
packing to come to a steady state temperature, the power 
supply was turned off and the multimeter was switched from 
DC voltmeter to ohmmeter. In this manner an accurate value 
of the paper resistance was determined at a temperature which 
corresponded to the liquid crystal event temperature. The 
multimeter was then switched back to DC voltmeter and the 
wind tunnel was started and brought up to speed. 

The wind tunnel velocity range and test section velocity 
profile were examined prior to obtaining experimental heat 
transfer data (see Appendix D). It was determined that the 
velocity profile was uniform at the test section entrance, 
and that the turbulence intensity was approximately 0.2%. 
After numerous checks with the hot wire anemometer, it was 
decided to use the installed wall manometer to set the tunnel 
speed of 98 ft/sec. 
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After the tunnel reached operating speed, the voltage 
was adjusted such that the liquid crystal color bands would 
appear at various angular positions around the cylinder. 

The first color band to appear would be at the lowest local 
Nusselt number on the cylinder. For tests with a d/r greater 
thpn 0.25 (where d = distance from plate to cylinder rim, and 
r = cylinder radius) this was always the laminar separation 
point (Figures 13-19). As voltage was increased, the color 
lines would appear in one, or as many as 5 or 6, angular 
locations all having equal Nusselt numbers. In this manner 
the cylinder surface temperature could be compared to the 
air stream temperature as determined by a thermocouple 
mounted in the wind tunnel stilling chamber. Using this 
temperature difference, and calculating the power dissipated 
from the voltage applied, the local heat transfer coefficient 
could be calculated for the circumference of the cylinder. 
Because of the way the cylinder was constructed, readings 
within 15 degrees of the seam were considered unreliable. 

The majority of Nusselt number data were collected with the 
seam located at 180 degrees. The Nusselt number in the 
vicinity of 180 degrees was determined by turning the cylinder 
such that the seam was facing forward at 0 degrees. Several 
experiments were also made with the seam at positions of 
135, 225, and 315 degrees. The final position of 0 degrees 
was chosen for the seam after it was observed that the value 
and angular location of the Nusselt number at the laminar 
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separation point, and in the vicinity of 130 and 210 degrees, 
was consistent with the values obtained with the seam at 
180 degrees. 

For all experiments, distinct, repeatable values of the 
Nusselt number versus angular location were obtained. The 
gradient as determined by the color bands, was sharp, with 
the exception of a region near the forward stagnation point, 
which appeared nearly isothermal. The gradient was also 
indistinct in the turbulent boundary layer, where the 
"scrubbing" action of cooler air caused the Nusselt number 
to change with time. Therefore, the Nusselt number shown 
for the turbulent regions must be considered to be an average 
value. 

B. PRESSURE COEFFICIENT EXPERIMENT 

With the pressure cylinder mounted in the wind tunnel, 
the manometer bank was leveled, pressure leads checked for 
tightness, and the wind tunnel brought up to a speed of 
98 ft/sec, corresponding to a Reynolds number of 153,000. 

The cylinder was positioned such that the pressure taps were 
initially at 0, 90, 180, and 270 degrees. The manometer 
readings for the four taps, the stilling chamber, and the 
test section were recorded. The cylinder was rotated in 
5 degree increments until pressure readings were obtained for 
the entire circumference. 
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The pressure coefficient was determined as follows: 





C P 


p-p p-p 

CO oo 
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whe re : 
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11 


pressure coefficient 




p 


density of air 




u 

oo 


air stream velocity (test section entrance) 




P 


pressure on cylinder surface 




P s = 


stagnation pressure (stilling chamber) 




p 

CO 


air stream static pressure (test section 
entrance ) 
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V. RESULTS 



Initial heat transfer and pressure coefficient experi- 
ments were made with the plate and cylinder attached 
(d/r = 0.0). A thin felt strip was placed between the two 
to assure a tight seal. The data obtained in this configura- 
tion agreed well with the results of McComas (see Figure 6). 
It can be seen that the trends in the data are identical, 
that local maximum and minimum values are within experimental 
uncertainty (see Appendix B for uncertainty analysis), and 
that a local maximum and local minimum were again observed 
near the point of attachment of the plate and cylinder. This 
is further evidence of the trapped vortex phenomenon as 
proposed by McComas. 

The next series of tests were made with the plate as 
far removed from the cylinder as possible, d/r = 5.33, or 
d = 8 inches. As there appeared to be some blockage, the 
method of Pope (4), as used by Meyer (5), was used to 
determine the necessary correction: 



Re 



corr 



(1.0 + e )Re 



oo 



where : 




and e 



sb 



is the solid blockage correction factor: 
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Figure 6. Comparison of the heat transfer results for d/r = 0.0 with 
the results of McComas at Re = 153,000. 




and e wb is the wake blockage correction factor: 

e wb * ’/4 <°/H) C D 

where : 

D = cylinder diameter, inches 
r = cylinder radius, inches 
H = total tunnel height, inches 
Cp = drag coefficient (see Meyer p. 107) 

Therefore : 

e - I (^l’ 2 + i <f?> d- 2 ) = °- 05 



and : 



Re corr = (1 - 05) 153 ’ 000 



Re 



corr 



161 ,000 



Using the R e corr to compare the results of this experi- 
ment with the work of Field, it was found that the results 
compare favorably with Field and with the theory of Schuh (6) 
(see Figure 7). It should be noted that Field experienced 
critical flow at Re = 149,000, while the experiments made in 
this work were exclusively subcritical for d/r = 5.33. This 
phenomenon could be attributed to higher free stream turbulence 
or cylinder deformation, either of which could have induced 
critical flow in Field's experiments. 
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Figure 7. Comparison of the heat transfer results 
for d/r = 5.33 and Re cQrr = 161,000 with 

the results of Field at Re = 148,000 
and with the theory of Schuh 
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Figures 13 through 19 are polar plots of the local 
Nusselt number and pressure coefficient for values of d/r 
ranging from 0.0 to 5.33. In examining the seven plots for 
significant trends or unusual phenomena, two items of 
interest were immediately apparent: one was the change from 
critical to subcritical flow, the other was a shift in the 
location of the stagnation point and the laminar separation 
points as the plate was moved away from the cylinder 
(see Figure 8) . 

The transition from critical to subcritical flow was 
investigated in order to ascertain the d/r ratio associated 
with the transition. As the flow was critical with the 
plate attached, and subcritical over the top at d/r = 0.25, 
the first transition point of interest was located in the 
interval 0.0 £ d/r <_ 0.25. By alternately moving the plate 
and adjusting the voltage, the transition point from critical 
to subcritical flow over the top of the cylinder was found 
to be approximately 0.08 £ d/r £ 0.10. 

As the flow over the bottom of the cylinder underwent 
a transition from critical to subcritical flow in the interval 
1.0 < d/r < 1.5, this region was investigated further. Using 
the same procedure of carefully adjusting spacing and voltage, 
the transition was found to occur at approximately 
1.17 < d/r < 1.25. 
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It should be noted that while both flows described 
above are definitely critical, there is significant 
difference in the degree of criticality on the top and 
bottom of the cylinder. This can be seen from the magnitudes 
of the Nu on the data shown in Figures 9 and 10. 

In order to glean further information from the data, 
the average Nusselt number, NlT, was computed from the 
equation: 

2 it 

Nu = ^ f Nu de 
0 

which was numerically integrated by computer using Simpson's 
Rule of Quadrature with a step size of 5 degrees. The 
results are shown in Figure 11a. 

It is interesting to note the unexpected shape of the 
curve. Nu" drops to a minimum value, then increases to a 
maximum, and then gradually decreases to a nearly constant 
value as d/r approaches relative infinity. These results 
prompted further investigation. 

In order to help determine the influence of the local 
Nu on Nu" a plot was made of Nu „ for the forward (-90 to 
90 degrees) and rear (90 to -90 degrees) portions of the 
cylinder. The results are shown in Figure lib and exhibit 
the same general shape as NTT. The drop in Nu from d/r = 0.0 
to 0.25 is explained from this plot when the sharp drop in 
N u m ax on the front of the cylinder from d/r = 0.0 to 0.25 is 
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d/r - 0.0 
d/r « 0.25 • 




Figure 9. Comparison of the local Nusselt number at 
d/r 0.0 with the local Nusselt number at 
d/r = 0.25. 
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Figure 10. Comparison of the pressure coefficient at 
d/r = 0.0 with the pressure coefficient at 
d/r = 0.25. 
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Figure lib. Maximum Nusselt numbers on the front and 
and back of the cylinder versus d/r. 
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noted. (The maximum local Nu for d/r = 0.0 occurs in the 
separation bubble, in the vicinity of 90 degrees, and is not 
used in the plot of Mu max -) 

Further investigation of the NU curve lead to an attempt 
to see how the voltage applied to the cylinder would vary if 
one location on the cylinder was maintained at a constant 
temperature as the plate spacing was varied. The locations 
150 to 165 degrees were chosen for this experiment. For each 
d/r the voltage was adjusted until a distinct color band 
appeared at these locations. Plotting the voltage applied 
versus d/r resulted in a curve of similar shape to NU and 

Nu max ^ see Fl ’9 ure 12). 

In order to verify the curvature of the Nu plot, a 
complete set of data was taken at d/r = 3.0. This data is 
not presented on a separate polar plot but NU was calculated 
and the value obtained enabled the shape of the Nu curve to 
be verified. 

A possible explanation for the observed variation in 
NU could be related to the amount of "cold" free stream air 
which has access to the cylinder surface. At d/r = 0.0, the 
flow is critical due to the fact that all of the air approach- 
ing the cylinder is forced to go "up and over." As soon as 
the plate is separated the flow has two paths. As soon as 
space is sufficient to provide a favorable path without large 
pressure gradients, enough of the air approaching the cylinder 
passes beneath to cause the flow over the top to transition 
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Vol tage vs d/r 
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Figure 12. Voltage required to maintain a predetermined temperature 
at the positions of 150 and 165 degrees versus d/r. 



to subcritical. Now the air passing below the cylinder is 
experiencing critical flow, possibly due to the nozzle 
effect produced by the rapidly decreasing flow area below the 
cylinder. The heat transfer in this region does not increase 
due to the fact that much air has already been heated passing 
over the front of the cylinder (i.e., the thermal and 
hydrodynamic boundary layers are a significant fraction of 
the gap size). 

As d/r is increased beyond 0.25, the larger opening 
allows more cold air, outside the boundary layer, to approach 
the bottom of the cylinder. This cold air carries off some 
heat by conduction and also mixes with the warm air in the 
separation bubble and turbulent boundary layer, resulting in 
increasing TTTT . 

The gradual decrease in Nu from, a high value at 
d/r = 2.0 is not understood at this time. It was expected 
that Nu would approach some constant value as the influence 
of the plate was reduced. It can be seen from Figure 11a 
that the uncertainty band nearly encompasses the decrease 
in Nu. Although the behavior of the Nu mgx and voltage curves 
suggest that the trend may be genuine, further investigation 
is needed to determine if Nu actually has a peak, or in fact 
approaches a constant value from some minimum point at a 
small d/r. 
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Figure 13. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 0.0. 
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Figure 14. 



Polar plot of Nusselt number 
coefficient versus angle for 



and pressure 
d / r = 0.25. 
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Figure 15. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 0.5. 



43 



091 



■ 350 . 



Nu # 
Cp * 




Figure 16. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 1.0. 
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Figure 17. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 2.0. 
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Figure 18 . 



Polar plot of Nusselt number and 
coefficient versus angle for d/r 



pressure 
= 4 . 0 . 



46 



180 



Nu © 



Cp x 




Figure 19. Polar plot of Nusselt number and pressure 
coefficient versus angle for d/r = 5.33. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The primary objective of this work, to design, construct 
and test an experimental apparatus which will allow the study 
of ground effects on a heated cylinder placed in a cross flow 
of air, has been achieved. The comparisons made with 
previous work are within experimental uncertainty, and 
indicate that the results of further experiments made with 
this apparatus can be relied upon. The liquid crystal 
technique for heat transfer analysis has proved to be an 
excellent method for determining quantitative heat transfer 
data. The results of the one plate geometry and Reynolds 
number investigated in this work are shown in the plots of 
Nu and Cp versus theta. Figure 8 shows the smooth shift 
in the location of key points which characterize flow 
phenomena (stagnation point, and laminar separation points) 
and makes it possible to approximate the value of Nu for 
configurations other than those investigated. 

Future work in this area should include a further 
investigation of the variation of NU with increasing values 
of d/r. The first step in this investigation could be 
further experi mentation with the configuration described in 
this report at various Reynolds numbers. Additional work 
should also investigate the effect of the plate length on 
flow patterns. As the wind tunnel used in this work has a 
capability for flow visualization, a series of experiments to 
determine flow patterns visually would prove to be valuable. 



48 



APPENDIX A 



HEAT TRANSFER AND PRESSURE DATA REDUCTION 

The equations and procedures used to reduce the heat 
transfer and pressure data are shown below. The results of 
each experiment are presented on polar plots of Nu and Cp 
versus Theta in section V. 



A. HEAT TRANSFER EXPERIMENT 

The relationship between power dissipated (by convection 
and radiation heat transfer), and power input by Joulean 
heating is the same as developed by Field (2) and used by 
Field and McComas ( 1 ) . 



h 



t 



RA(T-TJ (0.293) 



where: h. = total heat transfer coefficient, ~ 

1 hr ft ° F 

V = voltage applied, volts 
R = total resistance of Temsheet, Ohms 

2 

A = total heated surface area of Temsheet, Ft 

T = surface temperature of the cylinder as 
determined by liquid crystals, °F 

T w = air stream temperature, °F 



n watt-hr _ 746 watts 

BTU 2545 BTU/hr 
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The Nusselt number is defined as: 



Nu = 



h d 
c 



where : 



Nu 

h c 

d 

k 



= Nusselt number 



= convection heat transfer coefficient, 



= cylinder diameter, ft 
= 0.015 ^ r - ^Y<rp - (air at 75°) 



BTU 
hr ft 



2o, 



Since Nu is defined by the convective heat transfer 
coefficient, it is necessary to subtract the radiation effects 
from the total heat transfer coefficient 



This correction is made using the equivalent radiation heat 
transfer coefficient h^, where h^ is defined as: 

h = a(T + T ) (T 2 + T 2 ) 

y* \ — C oo oo 

where : 

h = equivalent radiation heat transfer 

r . * BTU 

coefficient * — 

hr f t ° F 

a = 0.1714 x 10" 8 ^ — x (Stef an-Bol tzman Const) 

hr ft' °R 

o 4 ,_2 = 0-9 radiation exchange factor between the 
cylinder and the surroundings 

T = cylinder surface temperature, °R 

T^ = air stream temperature, °R 
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A sample calculation for data obtained at a Re = 153,000 



and d/r = 0.50 is shown below: 



e = 68 degrees, V = 24.65 volts, T = 111.0°F, T ro = 63.5°F 
h. 



24. 65 2 volt 2 



t (13. 2*0 (0.314: f t 2 ) ( 1 1 1 -63 . 5° F ) ( . 293 ~ — | r g ^) 
BTU 



h t = 10.5 



hr ft °F 



hi = (0.1 714x1 O -8 ~ — x) ( • 9 ) ( 571 + 523. 5°R) 

r hr ft^°R 4 

( 571 2 + 523. 5 2 R 2 ) 



h r = 1.0 



BTU 



hr ft ° F 



V h r 



h _ = 9.5 



BTU 



hr f t 2 ° F 
BTU 



(9.5 — D - ' - u - 9 — )(.25 ft) 

I Y* A. £ o r- 



Nu = 



hr ft ° F 



(.° 15 ° hr B f to F ) 



Nu = 158 



B. PRESSURE COEFFICIENT EXPERIMENT 

As described previously, the pressure coefficient is 
defined as: 



C 



P 



P-P 



h p U‘ 



P-P 

oo 
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whe re : 



0^ = pressure coefficient 
p = density of air, lb/ft^ 

= air stream vel oci ty, ft/sec 
P = pressure on cylinder surface, in 



P $ = stagnation pressure, in 
P = air stream pressure, in 



h 2 o 



H„0 



h 2 o 



A sample calculation for data obtained at a Re = 153,000 
and d/r = 0.50 is shown below: 



e = 70° P = 10.0 P =6.0 P = 1 .86 

00 s 



10.1 - 6.0 
p 1.86 - 6.0 



C p = - 0.99 
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APPENDIX B 



UNCERTAINTY ANALYSIS 

Because of the uncertainties in the measured variables, 
it was necessary to perform an uncertainty analysis in 
order to ascertain the degree of uncertainty in the experi- 
mental results. The method of Kline and McClintock (7) 
was used. 

The following estimated uncertainties existed: 



voltage, V 


+ 0.1 


vol ts 


(20:1 ) 


resistance, 


R + 


0.2 ohms 


(20:1) 


temperature 


AT + 


_ 1 .0°F 


(20:1) 


area 


A + 


0.01 ft 2 


(20:1) 


diameter 


d + 


0.001 ft 


(20:1) 


radiation 
exchange factor 


S* + 


0.05 


(20:1 ) 


wall manometer 


H + 


0.05 in 


(20:1 ) 


manometer bank 


S + 


0.05 in 


(20:1 ) 



Total Heat Transfer Coefficient, h^: 

As shown in Appendix A, the equation used to calculate 
the total heat transfer coefficient was: 

h - - v 2 

h t " RA (T-T ) 

CO 

(since the conversion from watts to BTUs is a constant, 
it does not enter into the uncertainty analysis). 
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- 







For data obtained at a Reynolds number of 153,000 



and d/r = 0.50: 

V = 24.6/ + 0. 1 volts 
R = 13.2 + 0.3 ohms 
AT = 47.5+1.0 °F 
A = 0.315 + 0.001 ft 2 




— ~ =0.03 or to, = + 0.3 
h t h t " 

h t = 10.5 + 0.3 

In some regions where T = +_ 2.0 °F, then: 

to . 

jp fc = 0.05 or 0 ) h = 0.5 
h = 10.5 + 0.5 

All measured quantities were read with the uncertainties 
shown above, with the result that the calculated value of 
total heat transfer coefficient had a relatively constant 
value around the circumference of the cylinder. 
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Since the angular location of the color bands does not 
enter into the calculation of H^, the angular uncertainty 
has not been included here. This uncertainty is estimated 
to be + 3 degrees . 

Equivalent Radiation Heat Transfer Coefficient , h : 

The equation used to calculate was: 

li = ocR T + T )(T 2 + T 2 ) 

y* ' oo co 

defining T m = 1/2 (T + T m ) 
h r - 4o?TT l 

since a is a constant, it does not enter into the uncertain- 
ty analysis. 

For the same data used to perform the uncertainty 

analysis for h., the uncertainty in h was found to be: 
v r 




li = 1.00 + 0.06 
r — 
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Nusselt Number, Nu : 



The equation used to calculate Nu was: 
h d 

Nu = 

where h _ = h . - h 
c t r 

since k is a property of air, and its uncertainty is 
negl i gi bl e : 




Nu = 158+5 

“(In 

While this analysis has shown that = 0.03, and 
e = + 3 degrees, it should be noted that this may be 
misleading. Because there are portions of the cylinder which 
exhibit different flow phenomena, there are different types 
of color bands exhibited. The isothermal region has no 
distinct bands, and the mottled colors which appear cover 
a wide area. The liquid crystals were calibrated with 
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respect to the onset of each color. The result is that 
for areas of pulsating or mottled colors, the uncertainty 
of surface temperature and location are much higher then 
forregions in which a distinct color gradient exists. 



Reynolds Number , Re : 

The equation used to calculate Re was: 

U d 

_ co 

Re = 

v 

The value of is directly proportional to the sqaure 
root of the wall manometer reading, therefore: 



since v is a property of air, and its uncertainty is 
negligible: 



“Re 

Re 






2 



“Re 

Re 



0.011 or to Do = 1 700 *= 2000 
Re 



Re = 153,000 + 2000 
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Pressure Coefficient , Cp : 

The equation used to calculate C 



P 



was : 



C 



P 




A P 
TP 



1 

2 



both pressure differences were measured on the same manometer 
bank. For Re = 153,000 and d/r = 0.50: 



P = 10.1 P =6.0 P = 1 .86 

°° s 




C p =-0.99 + 0.02 
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APPENDIX C 



HOT WIRE ANEMOMETER CALIBRATION 

Prior to any experimental runs, the TSI hot wire 
anemometer was calibrated using a controlled air stream and 
a precision manometer. In order to assure accurate calibra- 
tion, the probe was mounted in the wind tunnel traversing 
mechanism and the calibration was done in the wind tunnel 
test section. In this way all leads and connections used 
for experimental runs would be in the circuit for calibration. 

The calibration apparatus provided a controlled air 
flow whose velocity was calculated from the pressure 
differential provided by an al cohol -f i 1 1 ed manometer. The 
velocity calculated was plotted against the voltage indicated 
by the hot wire anomometer unit (Figure 20). Calibration 
was done by increasing the calibration apparatus air stream 
velocity so that the alcohol in the manometer would indicate 
even values (in mm). The manometer inclination and readings 
were recorded, along with the voltage indicated by the hot 
wire anemometer. 

The equations used, data taken and calibration curves 
are shown below. Calibration checks were made periodically 
throughout the course of this work, and corrections made to 
account for the aging of the probe. 
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From Bernoulli's Equation: 



U 



■V 



2 9 c ap 



whe re : 

U = air stream velocity ft/sec 

g c = 32.17 ft/sec 2 

p = air density lb/ft 2 

2 

AP = pressure differential lb/ft 



The pressure differential from the manometer was calculated 
from the following relations: 



and 




g 



h 



0.789 

0.80 



ch ' 



where : 



hence : 



where : 



0.789 

0.80 



c 

h ' 



ratio of specific gravity of manometer 
fluid used, to specific gravity of 
fluid intended for manometer 

manometer inclination factor 

height of miniscus in inclined manometer mm 



U = 1 3. 24 /FT ft/sec 

co 



h = mm of alcohol 
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Turbulence intensity was determined by: 
-r 4 V e 



u 


v z -vg 




where: 






T u 


= turbulence intensity 




V 


= voltage indicated by hot wire 


anemometer 


e 


= rms voltage variation 




v o 


= extrapolated value of voltage 
to zero velocity 


correspond!' ng 
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APPENDIX D 



WIND TUNNEL VELOCITY PROFILE AND TURBULENCE INTENSITY 

In order to assure that test results would be accurate, 
it was necessary to show that the velocity profile in the 
test section was constant and that turbulence was small 
enough to have no effect on test results. 

The hot wire anemometer was used to probe the test 
section at various locations. The test section coordinate 
axis is shown in Figure 21. The test cylinder was located 
at X = 13.71, Y = -2.33. Velocity profiles were taken at 
the following locations: 

Test section entrance: X = 17.75 -9 . 0£Y^+9 . 0 Z = 0.0 

X = 17.75 Y = -2 . 2<_+4 . 0<_Z = 4.0 

Plate leading edge: X = -1.71 -3 . 023<_Y<_+9 . 0 Z = 0.0 

Top of the cylinder: X = -1 3.71 -0 . 825<_Y< + 9 . 0 Z = 0.0 

(plate attached) 

Top and bottom of cyl . : X = -13.71 -1 0 . 0<Y<_+9 . 0 Z = 0.0 

(plate at d/r = 5.33) 

All velocity profiles were taken at a tunnel speed of 
98 ft/sec. 

As shown in Figures 22 and 23 the air speed increased 
from 10 to 30 percent on either side of the cylinder in the 
case of the plate at d/r = 5.33, and a like amount over the 
top of the cylinder in the case of the attached plate. The 
velocity at the test section entrance was essentially constant. 
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Turbulence intensity was determined at the test section 



entrance: 



X = 17.7 - 5 . 0£ Y_<+5 . 0 Z = 0.0 

for three Reynolds numbers; 75,000, 100,000, and 153,000. 
The average value of turbulence intensity for Re = 153,000 
was 0.00239, or about 0.2%. The turbulence for the other 
Reynolds numbers was slightly higher, but of the same order 
of magnitude (see data). The value of turbulence intensity 
was calculated in the same manner as for the hot wire 
anemometer probe calibration described in Appendix C. 
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FLOW 

DIRECTION 



TEST SECTION TRANSLATOR 
pOORDINATE AXIS 

Figure 21. Isometric of the test section coordinate axis. 
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/ 7.7 0,0-1.11 -is.!/ 

Figure 22. Velocity profiles at the test section entrance, Plate 
leading edge, and the top of the cylinder for d/r = 0 
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Figure 23. Velocity profiles at the top and bottom of the cylindder for d/r 
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